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Abstract 

The transient response of an elastic cylinder to a laser impact is studied. When the laser 
source is a line perpendicular to the cylinder axis, modes guided along the cylinder are generated. 
For a millimetric steel cylinder up to ten narrow resonances can be locally detected by laser 
interferometry below 8 MHz. Most of these resonances correspond to Zero-Group Velocity guided 
modes while a few others can be ascribed to thickness modes. We observe that the theory describing 
the propagation of elastic waves in an isotropic cylinder is not sufficient to precisely predict the 
resonance spectrum. In fact, the texture of such elongated structure manifest as elastic anisotropy. 
Thus, a transverse isotropic (TI) model is used to calculate the dispersion curves and compare them 
with the measured one, obtained by moving the source along the cylinder. The five elastic constants 
of a TI cylinder are adjusted leading to a good agreement between measured and theoretical 
dispersion curves. Then, all the resonance frequencies are satisfactorily identified. 

PACS numbers: 43.20.Bi, 43.35.Cg, 43.20.Ks, 43.35.Zc 
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I. INTRODUCTION 


Elongated cylindrical structures like rods, cable strands or fibers are widely used in the 
industry or in civil engineering. Characterization of mechanical properties of constitutive 
materials is important for testing their structural integrity. Non Destructive Evaluation 
(NDE) of these properties is usually carried out with elastic waves. Various methods like 
Resonant Ultrasonic Spectroscopy (RUS)/pulse-echo contact-met hod, or pulsed/continuous 
laser-ultrasound (LU) contactless-method are used. Cylindrical structures support the prop¬ 
agation along their axes of elastic waves of different types: longitudinal (L), flexural (F) or 
torsional (T). Many theoretical and numerical studies and few experimental results were 
published this past fifty years on the propagation of time harmonic guided modes in solid 
or hollow cylinders. 1-3 In the nineteenth century, Pochhammer and Chree first established 
the equation of longitudinal waves in free isotropic cylindrical structures. In 1943, Hudson 
studied the fundamental flexural mode in a solid cylinder 4 and longitudinal modes of a bar 
were examined by Davies 5 in 1948. Gazis reported the first exact solutions of the frequency 
equation, 1 as well as a complete description of propagative modes, displacement and stress 
distributions for an isotropic elastic hollow cylinder in vacuum. Pao and Mindlin 6 ’ 7 as well 
as Onoe, McNiven and Mindlin 8 studied all the branches of the complete three dimensional 
problem of a free solid cylinder. A thorough review of elastic wave propagation in isotropic 
elastic cylinders and plates was given by Meeker and Meitzlcr. 3 Later, Zemanek investigated 
elastic wave propagation in a cylinder, both experimentally and theoretically. 9 

All these works deal with isotropic material, however transverse anisotropy is exhibited 
by elongated cylindrical structures due to their manufacturing processes (poly-crystalline 
metals) 10 or their texture (carbon fibers used in reinforced polymers). Then, non-destructive 
measurement of elastic constants of transversely isotropic (TI) materials is of great interest, 11 
especially in aeronautic and aerospace industries. Morse first established the frequency 
equation for longitudinal waves propagating in TI cylinders. 12 The extension of Gazis for¬ 
mulation to orthotropic and TI waveguides was initiated by Mirsky in 1964 for infinite and 
finite cylinders. 13 ’ 14 Then, several researchers investigated the propagation or the scattering 
of elastic waves in free or fluid-loaded TI cylinders. 15-19 Ahmad and Rahman 20 also studied 
the scattering of acoustic wave incident on a TI cylinder and showed that Buchwald’s repre¬ 
sentation 21 yields much simpler equations. 22 This representation is beneficial to simplify the 
description of the potential functions and economizes laborious calculations. This model 
provides perfectly similar results with Honarvar and Sinclair model 23,24 and can be applied 
to study both isotropic and TI cylinders. Honarvar et al. also obtained the frequency equa- 
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tions of axisymmetric and asymmetric free vibrations of finite TI cylinders. 25 Several other 
researchers proposed to use the impedance matrix theory 26,27 or a SAFE method to study 
transient thermoelastic waves in isotropic/anisotropic cylinders. 28 

Frequency equations determining the angular frequency oj(= 2irf) versus the axial wave 
number k(= 2tt/X ) have to be solved numerically. For a given circumferential order n, the 
various solutions (integer m) can be grouped into different families: longitudinal L(0, m) 
and torsional T(0,m) modes with displacements independent of the azimuthal angle 0, 
or flexural F (n,m) modes with displacements varying as sin n6 or cos n6. Each mode is 
represented by a dispersion curve oj(k). Many similarities exist between the families of 
modes found in plates and cylinders. Indeed, torsional modes T(0,m) are similar to shear 
modes SH in plates. Longitudinal L(0,m) and first flexural F(l,m) modes in cylinders are 
analogous to symmetrical S m and anti-symmetrical A m modes in plates. Major differences 
appear in the frequency spectrum for flexural modes F (n,m) with family number n ^ 2. 

Laser-based ultrasonic technique is a convenient tool for the generation and detection 
of guided elastic waves in plates 29 or cylinders. 30 It has been shown that this non-contact 
technique is very efficient to observe Zero Group Velocity (ZGV) Lamb modes corresponding 
to a frequency minimum of the dispersion curves oj(k). 31,32 For these specific points, the 
energy deposited by the laser pulse remains trapped under the source. The resulting local 
and narrow resonances can be detected at the epicenter with an optical interferometer. This 
non-contact method allowed us to perform accurate measurements of elastic properties of 
isotropic or transverse isotropic materials. 33,34 For cylinders, a line laser source can be used 
to control the propagation direction. For a line source parallel to the cylinder axis, cir¬ 
cumferential Rayleigh and whispering gallery waves are excited giving rise to resonances at 
all n integer values of the normalized circumferential wavenumber 35 ka. This configuration 
was also used by Mounier et al. to study resonances of a micrometric fiber at sub-gigahertz 
frequencies. 36,3 ^ 

In this paper, we investigate the mechanical response of a transversely isotropic cylinder 
to a laser line source perpendicular to the axis. In this case, axially guided elastic waves are 
generated and many resonances are observed. We performed measurements in a millimeter 
steel cylinder. Many resonances are observed in the MHz range. They can be roughly iden¬ 
tified from the minimum frequency of some branches of the dispersion curves of longitudinal 
and flexural modes calculated for an isotropic cylinder. However, significant discrepancies 
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remain between the isotropic model and experimental resonance frequencies. Dispersion 
curves are then calculated using a transverse isotropic model. The five independent elastic 
constants are adjusted such that all measured resonance frequencies can be precisely pre¬ 
dicted. Moreover, experimental dispersion curves, measured by moving the laser source, are 
compared with the theoretical ones. 


II. EXPERIMENTS 

Measurements were performed with an optical interferometer at the center of the source 
with the laser line perpendicular to the cylinder axis. Resonances were extracted from the 
spectrum of the temporal signal. Their frequencies were compared with minimum frequencies 
of dispersion curves calculated from an isotropic material. 


A. Laser-Ultrasonic setup 

The experimental setup is shown Fig. 1. The sample is an austenitic stainless steel (AISI 
304L) solid cylinder (length 420 mm, diameter 2 a = 0.775 mm). The longitudinal and 
transverse velocities are Vl = 5650 m/s and Vr = 3010 m/s respectively, and the measured 
mass density is p = 7.91 x 10 3 kg-m -3 . The rod is supported by two bevelled metallic 
pieces to reduce the mechanical contact (friction). Elastic waves were generated by a Q- 
switched Nd:YAG (yttrium aluminium garnet) laser providing 8-ns pulses of 15-rnJ energy 
at a 100-Hz repetition rate (Quantel Centurion). The spot diameter of the unfocused beam 
is equal to 2.5 mm. A beam expander (x4) and a cylindrical lens (focal length 250 mm) 
were used to enlarge and focus the laser beam into a narrow line on the surface. The optical 
energy distribution was close to a Gaussian and the absorbed power density was below the 
ablation threshold. The full length of the source at 1/e of the maximum value was found 
to be 10 mm and the width was estimated to be 0.3 mm. In the thermoelastic regime, the 
source is equivalent to a set of force dipoles distributed on the surface perpendicularly to 
the line. 

Local vibrations were measured by a heterodyne interferometer equipped with a 100-mW 
frequency doubled Nd:YAG laser (optical wavelength A= 532 nm). This interferometer is 
sensitive to any phase shift A0 along the path of the optical probe beam, and then to 
the mechanical displacement u r normal to the surface. The calibration factor (85 nm/V), 
deduced from the phase modulation A</> = 2iTu r /A of the reflected beam, was constant over 
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FIG. 1. Source and probe geometry used to excite and to detect local resonances in a cylinder. S: 
laser source, I: interferometer, C: cylindrical lens, (color online) 

the detection bandwidth (20 kHz - 20 MHz). Large low frequency phase-shifts due to thermal 
effect or first flexural F(l,l) mode are eliminated by interposing, before amplification, a high- 
pass filter having a cut-off frequency equal to 1.5 MHz. Measurements were conducted at 
room temperature (21 ±0.5°C). Signals detected by the optical probe were fed into a digital 
sampling oscilloscope and transferred to a computer. 


B. Zero-Group Velocity resonances 

The propagation of elastic waves along the cylinder axis is represented by dispersion 
curves calculated from Zemanek equation/’ The numerical algorithm used to find the roots 
of the secular equation was proposed by Seco and Jimenez. 38 Firstly, the cut-off frequencies 
are evaluated with a bisection method, then a zero-finding algorithm is applied to determine 
each branch successively. The roots were obtained with an acceptable error of less than 
10~ 6 . Finally, the dispersion curves are determined for longitudinal modes and the first seven 
families of flexural modes in less than one minute on a personal computer, with a resolution 
A k = 10 -3 mm -1 . As shown in Fig. 2(a), dispersion curves obtained for longitudinal L(0, m) 
and first flexural F(l,m) modes are similar to those obtained for an isotropic elastic plate. 
Dispersion curves of higher order flexural modes F(n, m) are plotted in Fig. 2(b). As the in¬ 
terferometer is only sensitive to the normal displacement, torsional modes are not presented. 

A typical signal, corresponding to the mechanical displacement normal to the cylinder 
surface, is given in Fig. 3(a). As previously explained, the oscillations in the first 5 ps are 
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L(0,m) and V(\,m) F(n,m) lor n = 2-6 



(a) (b) 


FIG. 2. Dispersion curves for an isotropic stainless steel solid cylinder (diameter 0.775 mm), a) 
Longitudinal L(0, m) and first flexural F(l,m) modes, b) Higher order flexural modes F(n,m). 
Minimum frequencies corresponding to ZGV modes are indicated by vertical arrows, (color online) 


due to the large displacements associated to the low frequency components of the flexural 
mode F(l, 1) similar to the Ao mode in plates. As shown in the insert, the low amplitude 
tail for t > 5 ps is not noise but high frequency oscillations due to the ZGV resonances. The 
spectra of out-off plane displacement are shown in Fig. 3(b). Seven resonances dominate 
between 0 and 8 MHz. The peak at 2.77 MHz can be ascribed to the ZGV resonance at the 
minimum frequency of the F(2, l)-mode, while peaks at 4.38 and 4.29 MHz correspond to 
L(0, 2) and F(3,1)-ZGV modes, respectively. The peak at 3.42 MHz can be associated with 
a thickness resonance [horizontal arrow in Fig. 2(a)] at the cut-off frequency (3.46 MHz) of 
the F(l, 3) mode. The peak at 5.61 MHz is relatively close to the minimum frequency (5.79 
MHz) of the F(4,1) mode. Higher frequency peaks at 6.83 and 7.70 MHz do not correspond 
to any minimum frequency on the dispersion curves. Except for L(0,2) which is estimated 
at 0.2%, relative errors from 2% to 4% are observed for the other resonances. 
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Thus, the isotropic model is not accurate enough: in the following section we apply a 
model developed by Ahmad et al. for transverse isotropic cylinder. 20 



(a) 



(b) 

FIG. 3. a) Radial displacement measured by the optical interferometer. Inset: zoom on a portion 
of the temporal displacement, b) Frequency spectrum filtered by an analog high-pass filter (/ ^ 1.5 
MHz), (color online) 

The Q-factor (Q = / 0 /A/ 0 ) can be estimated from the half-power width A/ 0 of the 
resonance peak at the ZGV point. In order not to underestimate the Q-factor, the signal 
acquisition time window 0 must be larger than the inverse of the bandwidth. We performed 
a measurement of the out-of-plane displacement at epicenter with a signal acquisition time 
window 0 = 4 ms. We obtained thinner resonances and a Q-factor at L(0,2)-ZGV frequency 
of 4 x 10 3 (A/ = 1.1 kHz and / 0 = 4.38 MHz). For flexural ZGV resonances, the Q-factor 
are slightly lower and vary from 2 x 10 3 to 2.5 x 10 3 . 
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III. TRANSVERSE ISOTROPIC MODEL 

In the linear theory of elasticity, anisotropic media are described by the stiffness tensor 
Cijki (with i, j, k, l = 1 to 3). Using the Voigt’s notation, they are represented by a 6 x 6 
symmetric matrix c af 3 (a, /3 — 1 to 6 ). Given the Cartesian coordinates (xi, x 2 , x 3 ) with 
the x 3 -axis parallel to the cylinder z-axis, the sample is supposed to be isotropic in the 
(xi, x 2 ) plane. Such transverse isotropic medium is described by five independent elastic 
constants : cn, C13 , C33 , C44 , and cqq. Other elastic constants are related to these coefficients: 
C22 = cn, c 23 = C13 , c 55 = c 44 , C12 = cn-2c 66 or vanish. 

A. Dispersion equation 

As mentioned by Honarvar et al. , 24 the displacement vector u{r, 0,z,t) can be derived 
from three scalar potential functions p, y, and -0. 23,39 A simple representation initiated by 
Buchwald 21 and used by Ahmad and Rahman 20 is the following, 



( 1 ) 


In cylindrical coordinates (r, 6 , z), the above representation leads to the following displace¬ 
ment components, 



{u r , Uq, U z ) 


Harmonic solutions are given by 


tp n = &J n (/3r) cos (nd) exp [i(kz — cot)], 
0 ra = d' J n ((3r) cos (nd) exp [i{kz — cot)], 
Xn = XJ n (/3r) sm(n9) exp [i{kz — cot)], 


where J n is the Bessel function of first kind of order n. They correspond to the super¬ 
position of plane waves of wavevector k = [j3 cos 6, ft sin 6, k) which satisfy the propagation 
equation in a meridian plane. The radial component /3 of the wavevector k must satisfy the 
Christoffel equation 


(C 11 C 44 / 3 4 — A 0 2 + F ) (c66/3 2 + C 44 A : 2 — pco 2 ) — 0 . 


( 2 ) 


The first term correspond to waves polarized in the meridian plane (x 2 , x 3 ) and the second 
term to pure shear wave (xi). The coefficients E and F are given in Appendix [Eq. (A.2)]. 


The above equations are similar to those obtained by Mirsky. 14 Omitting the propagation 
term exp i(kz — cot) for simplicity, the three independent solutions are found to be 

(<pin, Vhn, xin) = ^J n (M cos(n9), q 1 J n (ftr) cos(n6>), o), 

(<^2 n, 4>2m X 2 n) = [q .2 Jnifor) cos(n6>), J n (for) cos(nd ), 0) , 

((P3n,^3n,X3n) = (o, 0, J n {far) sin(n6>) j , 

where q± and q 2 are the potential amplitude ratios provided in Appendix Eq. (A. 4). Hence, 
the displacements are as follows 

(uir, Uw, u lz ) n = (fiiJ' n (for) cos(n0), - J n ((3ir)sm(n9), ikq x J n (for) cos(n0)) , 

(U 2 r, “ 20 , u 2z )n = (q 2 foJ' n (for) cos(nQ), - 0) q 2 J n (for) sin(n6>), ikJ n (for) cos(nd) j ,(3) 

(« 3 r, “30, “3*)n = ( 0 ) J n (M cos(nd), -/ 3 3 J' n (for)sin(n6 ), o) . 

The first two solutions correspond to a coupling between quasi-longitudinal and quasi-shear 
waves polarized in the meridian plane. The last solution having zero displacement component 
along the 2 -axis corresponds to a pure shear wave. These solutions must be combined with 
weighting factors B n , C n , D n to satisfy the boundary conditions at the free surface r = a. 
These conditions imply that normal stresses vanish: a rr (a ) = a r g(a) = a rz (a) = 0 and lead 
to the following homogeneous linear system 


On Oi2 Oi3 


B n 


0 

<221 «22 «23 


c n 

= 

0 

O 31 O32 o 33 


D n 


0 


where the matrix elements are listed in the Appendix Eq. (A.l). The dispersion equation 
results from the secular equation 

det(dij) = 0. (5) 

For n — 0, Eq. (5) splits into two parts: a 23 = 0, he., (foa) Jq (foa) = J' 0 (foa ) similar to the 
equation given by Mirsky 14 for torsional modes T(0,m) and ana 32 = a 32 a 3 i corresponding 
to longitudinal modes L(0,m). 

B. Cutoff frequencies 

As the wavenumber k is equal to zero, the motion at the cutoff frequencies is independent 
of the axial coordinate x 3 . This implies k 2 qi = 0, q 2 = 0 and a\ 2 = a 22 = 0. The dispersion 
equation [Eq. (5)] for longitudinal and flexural modes in a solid cylinder at k = 0 reduces to 

fooJ n (foa)(a\ 3 a 2 i — ai±a 23 ) = 0, (6) 


9 








where the four a tJ simplify as 


an = cn(/3ia) 2 J"(/3ia) + ci 2 {(/?ia) J' n ((3ia) - n 2 J n (f3ia )}, 

°i 3 = 2 nc m {(/3 3 a)J' n (/3 3 a) - J n (/3 3 a)}, 

021 = -2n{(/3ia)4(^ia) - J n (£ia)}, 

a 2 3 = -(/3 3 a) 2 J"(/3 3 a) + (fi 3 a)J' n (fi 3 a) - n 2 J n (/3 3 a). 

To find the cutoff frequencies, the above equation can be solved for cu c (= 27r/ c ), for various 
values of rt and m. The hrst two roots /?i 2 are simplihed using the reduced expressions of 
A = (cn — c 44 )pu; 2 , E = (cn + c 44 )/xu 2 and F = p 2 u^ [Eq. (A.3)] as /3 3 = oj c ^p/cu and 
@2 = u c \Jp/c 44 . The third root reduces to f3 3 = u c \Jp/c§§. The cutoff frequency equation 
[Eq. (6)] is simplified as = 0 for L(0,m). 


C. Dependence of ZGV frequencies on elastic constants 

A sensitivity analysis was performed to characterize the elastic constant influence on the 
dispersion curves and especially on ZGV frequencies. Figure 4 displayed the longitudinal 
and the hrst two flexural mode families as a function of elastic constants. Each constant 
was successively varied by ±5% and ±10% around the following average values (in GPa): 
cn = 240, C 13 = 110, C 33 = 250, c 44 = 70, and cqq = 70. The curves obtained for the average 
values are displayed in red dashed line. It appears that each ZGV resonance depends on 
specific constants. For example, the elastic constants which mainly affect the L(0,2)-ZGV 
frequency are Cn and c 44 . These elastic constants are proportional to the square of lon¬ 
gitudinal and shear velocities in directions perpendicular to the cylinder axis x 3 . While 
F(1,4)-ZGV frequency slightly depends on Cn and more strongly on c 44 and c 66 which is 
proportional to the shear velocity in the direction perpendicular to the cylinder axis. When 
the values of Cn, c 43 and c 33 increase or when c 44 , cqq decrease, the F(l,4) minimum ZGV 
frequency disappears, i.e., the slope of the spectral line at the cutoff frequency becomes pos¬ 
itive. Then, F(2,1)-ZGV frequency is almost independent of the hrst four elastic constants 
while it significantly depends on c 66 . 


IV. RESULTS AND DISCUSSION 

In this section, we identify all the resonance frequencies and accurately determine the 
Eve elastic constants. The theoretical dispersion curves obtained with the TI model 20 were 
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Wavenumber k (mm 1 ) 


FIG. 4. Influence of the elastic constants on the dispersion curves of a transversely isotropic 
cylinder. From top to bottom: for longitudinal modes and the first two flexural mode families. 
Red arrows depict the ZGV frequencies and black arrows show increasing c t j. (color online) 

fitted with those obtained experimentally. Then, we compared isotropic and TI models to 
highlight the anisotropy of the stainless steel cylinder. 


A. Dispersion curves measurements 

Experimental dispersion curves were measured with the laser ultrasound setup displayed 
in Fig. 1. The out-of-plane displacement was recorded with the interferometer, at a single 
point in the middle of the cylinder, while the line source was moved along the cylinder axis 
on 40 mm by 0.1 mm steps. For each source position, the normal displacement was recorded 
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during 180 ps at a 50 MHz sampling frequency. 1024 signals were averaged to increase the 
signal to noise ratio. First, apodization (Hanning) windows were applied in both dimensions 
(time and distance) to avoid secondary lobes. Then, a 2D-Fourier transform was applied to 
the obtained B-scans. In order to observe the backward modes, the spatial Fourier transform 
was calculated for negative and positive wavenumbers (k) as shown in Fig. 5(b). 
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FIG. 5. a) Spatial Fourier transform of the normal displacement at 4.40 MHz illustrating backward 
propagation with a negative wavenumber. Each backward mode is coupled to the upper mode of 
the same family (n,m + 1). To avoid confusion they are indicated by the index b. b) Experimental 
dispersion curves obtained for negative and positive wavenumbers, c) Negative dispersion curves 
Lo(—k) were folded and added to the positive ones. The experimental dispersion curves were fitted 
with the TI model (the control parameters used are listed in the Tab. I). (color online) 


As previously discussed for Lamb modes in plates, 32 


cylindrical waveguides support 
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backward-propagation with opposite group and phase velocities. 40 The counter propagative 
modes are also observed in the cylinder both for longitudinal and flexural modes and occur 
in the vicinity of k — 0. For example, at the L(0,2)-ZGV frequency (4.34 MHz) the L(0,2) 
and L(0,3) modes having opposite wave vectors interfere. The power spectrum was com¬ 
puted at a frequency of 4.40 MHz, roughly higher than the L(0,2)-ZGV resonance frequency, 
for which the modes are propagative. Figure 5(a), shows that this spectrum is composed 
of two main peaks. The peak at a negative value of k, similar to the larger one in the 
positive k domain demonstrates the backward propagation in cylinders. Afterwards, the 
dispersion curves (backward region) for negative wavenumber u>(—k) was folded and added 
to positive ones [Fig. 5(c)]. We can clearly identified five ZGV frequencies. One at 4.34 MHz 
corresponding to the L(0,2)-ZGV resonance and four others at 2.77, 4.28, 5.60 and 6.60 MHz 
corresponding to F(2,l), F(3,l), F(4,l) and F(l,4) flexural-ZGV resonances, respectively. 
The wavelength at these ZGV frequencies extend from 2.3 mm to 3.2 mm (X/d ~ 3 — 4) 
except for F(1,4)-ZGV frequency which has a 14.3-mm wavelength (X/d ~ 18.6). 


TABLE I. Elastic constants Cij (or sound velocities) used in the TI model, a and p were measured. 


Stiffness constants (GPa) 

cn 

C12 

C13 

C33 

C44 

cm 

239 

108 

110 

252 

72 

65 


Sound velocities ( 

m/s) 


Diameter 






(pm) 

Fli 

Vsi 

Els 

F53 

v b 

2 a 

5500 

2875 

5650 

3010 

4808 

775 ±2 

Young modulus (GPa) 


Poisson’s ratio 


Density 






(kg/m 3 ) 

E± 

E n 

V\2 

PL 3 

Z'.Sl 

P 

183 

172 

0.317 

0.316 

0.298 

7910 ± 5 


The theoretical dispersion relation u(k) [Eq. (5)] were calculated by solving the secular 
equation with the zero-finding algorithm. 38 Experimental dispersion curves were fitted with 
the Tl model. Bulk modes proportional to stiffness constants C33 and C44 were eliminated 
by dividing the determinant by (pu 2 — c 33 k 2 )(pu 2 — c 44 fc 2 ). We calculated longitudinal 
modes and the first seven families of flexural modes. The theoretical dispersion curves were 
fitted heuristically with five control parameters and with the diameter and the mass density 
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previously measured. We can observed [Fig. 5 (c)] a good agreement between experimental 
dispersion curves and those predicted by the TI model. The parameters used for the TI 
model are listed in Tab. I. The red and blue dashed lines correspond to the longitudinal and 
first flexural modes respectively, the green dashed lines indicated the flexural modes with 
n ^ 2. Although most modes are well fitted by the theory, a few experimental branches are 
not explained. This may be due to a small remaining error on the estimated constants or 
to a slight discrepancy with transverse anisotropy. 


B. Elastic parameters 


We now explain the physical meaning of the five sound velocities implicitly involved in 
the TI model. Stiffness constants Cn and C33 determine the longitudinal (L) ultrasound 
velocities in perpendicular (Vta) and parallel (Vl 3 ) directions with respect to the cylinder 
axis £3 

cn = pV'li , C33 = pV‘L- ( 7 ) 

Stiffness constants C44 and C66 = (cn — ci2)/2 determine the shear (S) velocities in directions 
perpendicular to the cylinder axis with polarizations either parallel (Vgi) or perpendicular 
(V S3 ) to x 3 -axis 

c 44 = pVfp c 66 = pVf 3 . (8) 


The remaining stiffness constant C13 should be adjusted manually. Hence, two Young’s 
modulus can be defined: E// — E 33 for a stress parallel to the cylinder axis and E± = En = 
E22 for a stress perpendicular to the cylinder axis. They are given by 


E// = 


Cll — Cqq 


-,e ± = 


4 c 66 c 2 


C + C33C66 


( 9 ) 


where c 2 = 033(044 — cqq) — cf 3 is the determinant of the sub-matrix of elastic constants c t j. 
This effective stiffness constant c is proportional to the bar velocity of the waves propagating 
along the cylinder axis: Vj, = {E///p ) 1 ^ 2 (see Mirsky 14 Eq. 44 ). Three Poisson’s ratios 41,42 
can be defined, 1/42 and z/ 31 for longitudinal extension in the basal plane and z/ 13 along the 
6-fold axis 

C42C33 — C 2 3 _ 2C66C43 _ C43 

v \2 — —5—;-, ^31 — -) ^13 — -;-• 

C Z + C33C66 C z + C33C66 C44 + C42 


C. Comparison between isotropic and TI model 

Theoretical dispersion curves obtained with isotropic (blue dashed lines) and transversely 
isotropic (red lines) models are plotted in Fig. 6. The longitudinal modes are displayed in 
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[Fig. 6(a)] and the first six flexural mode families [Figures 6(b)-6(g)]. 



FIG. 6. Comparison of isotropic model (blue dashed lines) with TI model (red lines), a) Lon¬ 
gitudinal mode L(0, m). b-g) Flexural modes F(n, m) for n extend from 1 to 6. h) Frequency 
spectrum recorded at epicenter. Inset: two resonance frequencies are observed corresponding to 
F(1,4)-ZGV and F c (4,2), respectively. The red square (□) indicate the cutoff frequencies calculated 
with [Eq. (6)]. (color online) 

The frequency spectrum recorded at the epicenter [Fig. 6(h)]. This representation allows 
to distinguish the different propagation modes and to precisely identify each resonance peak. 
It obviously appears that the isotropic model is not appropriate to interpret experimental 
results and clearly confirms the transverse isotropic nature of the cylinder. Five ZGV reso¬ 
nances, F(2,l), F(3,l), L(0,2), F(4,l), F(l,4), and eight thickness resonances, F(l,3), L(0,3), 
F(4,2), F(5,l), F(l,5), F(6,l), F(l,6) and F(2,5), can be identified from the dispersion curves 
calculated with the TI model. Both ZGV and cutoff frequencies are estimated with a relative 
error less than 0.6%, except for F c (l,3) which was estimated with an error close to 1.1%. 


D. Displacements at cutoff frequencies 

In order to explain why several thickness resonances (at k = 0) are observed while others 
are not, we analysed the displacements generated inside the cylinder. To this end, we 
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calculated u r , Uq and u z , using Eq. (3). The radial and axial displacements ( u r ) and ( u z ) 
were calculated at 6 = 0 and the azimuthal displacement ( uq ) at 6 max = n/2n. They are 
displayed in Fig. 7 at the cutoff frequencies of the modes L(0,3), L(0,4), F(l,2) and F(l,3) 
which respectively correspond to the frequencies 4.834, 8.673, 2.276 and 3.343 MHz. For 
L c (0,3) and F c (l,3) thickness resonances, the radial displacement at the cylinder surface 
r = a is significant. Conversely, for the modes L c (0,4) and F c (l,2) the radial displacement 
vanishes for r = a, but have essentially a large displacement along z-axis. This explains 
why, at these cutoff frequencies, no resonances are observed. 



FIG. 7. Azimuthal displacement produced at a) L(0,3), b) F(l,2), c) L(0,4) and d) F(l,3) cutoff 

frequencies. Black (—), red dash (-) and blue (-) dash lines indicate u r , uq and u z , respectively. 

(color online) 


V. CONCLUSION 

In summary, elastic guided waves propagating in a stainless steel cylinder of millimetric 
diameter were investigated by laser ultrasonic techniques. Using a laser line source per¬ 
pendicular to the cylinder axis, ZGV resonances were observed both for longitudinal and 
flexural modes. One of the major interests is to identify each resonance to characterize 
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clastic properties of materials. We demonstrated that isotropic model was not relevant to 
describe the behaviour of elastic guided waves in a steel rod. A transverse isotropic model 
allowed us to describe the propagation and to calculate dispersion curves co{k). The five 
elastic constants were estimated from experimental dispersion curves measured by the LU 
technique. A good agreement between theoretical and experimental dispersion curves was 
obtained both for longitudinal and flexural modes. 

In the future, it would be useful to build a numerical inversion procedure to improve 
the accuracy of the material parameter estimation. Furthermore, it will be interesting to 
investigate the relationship between ZGV modes and elastic constants. This could be done 
by calculating the second order derivatives of the dispersion curves at cutoff frequencies to 
determine the existence of backward modes. Finally, it will be interesting to conduct other 
studies in materials with different types of anisotropy. This technique can be also applied in 
the sub-gigahertz range with a laser line source of width comparable to the cylinder diameter 
to study micrometric fibers. 
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Appendix 


The scalar functions p, x and ip satisfy the following wave motion equations 


Cll (V 

(C13 + C44) 



d 2 ip , , d 2 ip 

v + C44 ^W + ( Cl3 + C44 ) ~QJ2 


d z 2 


0)6 


^ + C44 





d 2 ip 

+C33 a^ 

d 2 x 
X + 


P 


d V 


p 

p 


d 2 ip 

d 2 X 
dt 2 ' 


The longitudinal L (represented by ip) and the vertically polarized quasi-transverse SV 
(represented by ip) waves are coupled, whereas the pure transverse wave SH (represented 
by x) is decoupled from the others. The dispersion equation of torsional, longitudinal and 
flexural guided modes in a transversely isotropic cylinder with free boundary conditions 
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results in the vanishing of the determinant of a 3 x 3 matrix a^- 

a 23 = 0 for T(0,m), 

ana 32 - ai 2 a 3 i = 0 forL(0,m), 
det(ciij) = < an(a 22 a 33 — a 23 a 32 ) — 
ai 2 (a 2 ia 33 — a 23 a 3 i) + 

k ai 3 (a 2 ia 32 - a 22 a 3i ) = 0 for F (n > 1 ,m), 
where the matrix elements are 

flu = Cu(PiaY J n (P\a) + 

ci 2 |(/3ia)X(/3ia) - n 2 J n {(3 ia) j - 
Ci 3 gi(fca) 2 J n (/3ia), 
oi 2 = cnq 2 ((3 2 a)~ J n ((3 2 a) + 

ci 2 g 2 1 (f3 2 a)j' n (/3 2 a) - n 2 J n (/3 2 a) j - 
c 13 (ka) 2 J n (/3 2 a), 

a i3 = 2nc 66 | ((3 3 a)j' n ((3 3 a) - J n (/3 3 a)| , 

a 2i = —2n | (/?ia) J n {f3\a) - J n {/3 ia)|, (A.l) 

a 22 = -2 nq 2 { ((3 2 a) f n (f3 2 a) - J n (/3 2 a) j , 

a 23 = — {Y 3 aY J n ((3 3 a) + (/ 3 3 a)J n (ft 3 a ) — n 2 J n ((3 3 a), 

a 3 i = (1 + gi)(/3ia)X(/3ia), 

° 3 2 = (1 + Q2)(fd2<yJ n (fd 2 a), 
a 33 = nJ n (/3 3 a), 

respectively. The first and second derivatives of the Bessel function of the first kind of order 
n can be expressed in terms of J n _i, J n , J n + 1 etc. by using recurrence relations. 

The coefficients E and F used in the Christoffel Equation [Eq. 2] are defined as 

E = (ci 3 + c±f) 2 k 2 + C 44 (pen 2 — C 44 A; 2 ) + cn(pur — c 33 k 2 ) 

= A 2 + C 44 F 1 + cnF 2 , (A. 2 ) 

F = (pen 2 - c u k 2 )(pu 2 - c 33 k 2 ) = F 1 F 2 . 

Finally, the three roots (/3) of Christoffel equation are 


Pi, 2 — 


and P3=J—, 
2C11C44 V Cqq 


where A = y/E 2 — Acuc^F. At last, the amplitude ratio q\ and q 2 are given by 

F\ - CnPl kE 1 


Qi = 


kEi 


Q2 = 


F\ - c u P 2 


(A.3) 


(A.4) 


18 









1 D. C. Gazis, Three-dimensional investigation of the propagation of waves in hollow circular cylin¬ 
ders. I. Analytical foundation. II. Numerical results, J. Acoust. Soc. Am. 31, 568-578 (1959). 

2 R. Mindlin and H. McNiven, Axially symmetric waves in elastic rods, J. Appl. Mech. 27, 145— 
151 (1960). 

3 T. Meeker and A. Meitzler, Guided wave propagation in elongated cylinders and plates, 
Phys. Acoust. 1, 111-167 (1964). 

4 G. E. Hudson, Dispersion of elastic waves in solid circular cylinders, Phys. Rev. 63, 46 (1943). 

5 R. Davies, A critical study of the hopkinson pressure bar, Phil. Trans. Roy. Soc. London. A. 
pages 375-457 (1948). 

6 Y.-H. Pao and R. Mindlin, Dispersion of flexural waves in an elastic, circular cylinder, 
J. Appl. Mech. 27, 513-520 (1960). 

' Y.-H. Pao, The dispersion of flexural waves in an elastic circular cylinder, part 2, J. Appl. Mech. 
29, 61-64 (1962). 

8 M. Onoe, H. McNiven, and R. Mindlin, Dispersion of axially symmetric waves in elastic rods, 
J. Appl. Mech. 29, 729-734 (1962). 

9 J. Zemanek Jr, An experimental and theoretical investigation of elastic wave propagation in a 
cylinder, J. Acoust. Soc. Am. 51, 265-283 (1972). 

10 T. Mason, Variation in the dispersion of axisymmetric waves in infinite circular rods with 
crystallographic wire texture, J. Acoust. Soc. Am. 106, 1262-1270 (1999). 

11 R. G. Payton, Elastic wave propagation in transversely isotropic media, volume 4, Springer 
Science & Business Media (1983). 

12 R. Morse, Compressional waves along an anisotropic circular cylinder having hexagonal sym¬ 
metry, J. Acoust. Soc. Am. 26, 1018-1021 (1954). 

13 I. Mirsky, Axisymmetry vibrations of orthotropic cylinders, J. Acoust. Soc. Am. 36, 2106-2112 
(1964). 

14 I. Mirsky, Wave propagation in transversely isotropic circular cylinders, part I: Theory, part II: 
Numerical results, J. Acoust. Soc. Am. 37, 1016-1026 (1965). 

15 F. C. Eliot and G. Mott, Elastic waves propagation in circular cylinders having hexagonal 
crystal symmetry, J. Acoust. Soc. Am. 44, 423-430 (1968). 

16 H. McNiven and Y. Mengi, Dispersion of waves in transversely isotropic rods, 
J. Acoust. Soc. Am. 49, 229-236 (1971). 

17 P. B. Nagy, Longitudinal guided wave propagation in a transversely isotropic rod immersed in 


19 



fluid, J. Acoust. Soc. Am. 98, 454-457 (1995). 

18 M. J. Berliner and R. Solecki, Wave propagation in fluid-loaded, transversely isotropic cylinders, 
part I. Analytical formulation, part II. Numerical results, J. Acoust. Soc. Am. 99, 1841 1853 
(1996). 

19 A. J. Niklasson and S. K. Datta, Scattering by an infinite transversely isotropic cylinder in a 
transversely isotropic medium, Wave Motion 27, 169-185 (1998). 

20 F. Ahmad and A. Rahman, Acoustic scattering by transversely isotropic cylinders, 
Int. J. Eng. Sci. 38, 325-335 (2000). 

21 V. Buchwald, Rayleigh waves in transversely isotropic media, Q. J. Mech. Appl. Math. 14, 293- 
317 (1961). 

22 F. Ahmad, Guided waves in a transversely isotropic cylinder immersed in a fluid, 
J. Acoust. Soc. Am. 109, 886-890 (2001). 

23 F. Honarvar and A. Sinclair, Acoustic wave scattering from transversely isotropic cylinders, 
J. Acoust. Soc. Am. 100, 57-63 (1996). 

24 F. Honarvar, E. Enjilela, A. N. Sinclair, and S. A. Mirnezami, Wave propagation in transversely 
isotropic cylinders, Int. J. Sol. Struct. 44, 5236-5246 (2007). 

25 F. Honarvar, E. Enjilela, and A. N. Sinclair, Asymmetric and axisymmetric vibrations of finite 
transversely isotropic circular cylinders, Acoust. Phys. 55, 708-714 (2009). 

26 A. N. Norris and A. Shuvalov, Wave impedance matrices for cylindrically anisotropic radially 
inhomogeneous elastic solids, Q. J. Mech. Appli. Math. 63, 401-435 (2010). 

27 A. N. Norris, A. J. Nagy, and F. A. Amirkulova, Stable methods to solve the impedance matrix 
for radially inhomogeneous cylindrically anisotropic structures, J. Sound Vib. 332, 2520-2531 
(2013). 

28 R. Chitikireddy, S. K. Datta, A. H. Shah, and H. Bai, Transient thermoelastic waves in an 
anisotropic hollow cylinder due to localized heating, Int. J. Solids Struct. 48, 3063-3074 (2011). 

29 D. A. Hutchins, K. Lundgren, and S. B. Palmer, A laser study of transient Lamb waves in thin 
materials, J. Acoust. Soc. Am. 85, 1441 (1989). 

30 D. Gsell and J. Dual, Non-destructive evaluation of elastic material properties in anisotropic 
circular cylindrical structures, Ultrasonics 43, 123-132 (2004). 

31 C. Prada, O. Balogun, and T. Murray, Laser based ultrasonic generation and detection of 
Zero-Group Velocity Lamb waves in thin plates, Appl. Phys. Lett. 87, 194109 (2005). 

32 C. Prada, D. Clorennec, and D. Royer, Local vibration of an elastic plate and Zero-Group 
Velocity Lamb modes, J. Acoust. Soc. Am. 124, 203 (2008). 


20 


33 


D. Clorennec, C. Prada, and D. Royer, Local and noncontact measurements of bulk acous¬ 
tic wave velocities in thin isotropic plates and shells using zero group velocity Lamb modes, 
J. Appl. Phys. 101, 034908 (2007). 

34 M. Ces, D. Royer, and C. Prada, Characterization of mechanical properties of a hollow cylinder 
with zero group velocity Lamb modes, J. Acoust. Soc. Am. 132, 180 (2012). 

35 D. Clorennec, D. Royer, and H. Walaszek, Nondestructive evaluation of cylindrical parts using 
laser ultrasonics, Ultrasonics 40, 783-789 (2002). 

36 D. Mounier, C. Poilane, H. Khelfa, and P. Picart, Sub-gigahertz laser resonant ultrasound 
spectroscopy for the evaluation of elastic properties of micrometric fibers, Ultrasonics 54, 259- 
267 (2014). 

37 H. Khelfa, D. Mounier, C. Poilane, and P. Picart, Evidence of guided acoustic waves propagating 
along a micrometric fiber, Appl. Phys. Lett. 105, 161906 (2014). 

38 F. Seco and A. R. Jimenez, Modelling the generation and propagation of ultrasonic signals in 
cylindrical waveguides, Ultrasonic Waves, Dr Santos (Ed.), ISBN: 978-953-51-0201-4, InTech 
( 2012 ). 

39 Y. Pan, M. Perton, C. Rossignol, and B. Audoin, The transient response of a transversely 
isotropic cylinder under a laser point source impact, Ultrasonics 44, e823-e827 (2006). 

40 H. Cui, W. Lin, H. Zhang, X. Wang, and J. Trevelyan, Characteristics of group velocities of 
backward waves in a hollow cylinder, J. Acoust. Soc. Am. 135, 3398-3408 (2014). 

41 S. Tokmakova, Anisotropy of Poisson’s ratio in transversely isotropic rocks, in AIP Conf. Proc., 
volume 1022, pages 413-416 (2008). 

42 A. Ballato, Poisson’s ratios of auxetic and other technological materials, 
IEEE Trans. Ultra. Ferro. Freq. Cont. 57, 7-15 (2010). 


21 


LIST OF FIGURES 


1 Source and probe geometry used to excite and to detect local resonances in 

a cylinder. S: laser source, I: interferometer, C: cylindrical lens, (color online) 5 

2 Dispersion curves for an isotropic stainless steel solid cylinder (diameter 0.775 

mm), a) Longitudinal L(0,m) and first flexural F(l,m) modes, b) Higher 
order flexural modes F(n, m). Minimum frequencies corresponding to ZGV 
modes are indicated by vertical arrows, (color online) . 6 

3 a) Radial displacement measured by the optical interferometer. Inset: zoom 
on a portion of the temporal displacement, b) Frequency spectrum filtered 

by an analog high-pass filter (/ ^ 1.5 MHz), (color online). 7 

4 Influence of the elastic constants on the dispersion curves of a transversely 

isotropic cylinder. From top to bottom: for longitudinal modes and the first 
two flexural mode families. Red arrows depict the ZGV frequencies and black 
arrows show increasing c % j. (color online) . 11 

5 a) Spatial Fourier transform of the normal displacement at 4.40 MHz illustrat¬ 
ing backward propagation with a negative wavenumber. Each backward mode 
is coupled to the upper mode of the same family (n,m+l). To avoid confusion 
they are indicated by the index b. b) Experimental dispersion curves obtained 
for negative and positive wavenumbers, c) Negative dispersion curves uj(—k) 
were folded and added to the positive ones. The experimental dispersion 
curves were fitted with the TI model (the control parameters used are listed 

in the Tab. I). (color online) . 12 

6 Comparison of isotropic model (blue dashed lines) with TI model (red lines), 
a) Longitudinal mode L(0,m). b-g) Flexural modes F(n, m) for n extend 
from 1 to 6. h) Frequency spectrum recorded at epicenter. Inset: two res¬ 
onance frequencies are observed corresponding to F(1,4)-ZGV and F c (4,2), 
respectively. The red square (□) indicate the cutoff frequencies calculated 


with [Eq. (6)]. (color online) . 15 

7 Azimuthal displacement produced at a) L(0,3), b) F(l,2), c) L(0,4) and d) 

F(l,3) cutoff frequencies. Black (—), red dash ( - ) and blue ( - ) dash 

lines indicate u ri uq and u z , respectively, (color online) . 16 


22 








LIST OF TABLES 


I Elastic constants c l3 (or sound velocities) used in the TI model, a and p were 
measured. 


13 


23 



